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The hole-state random phase approximation (hRPA) and the particle-state random phase approx- 
imation (pRPA) for systems like odd A nuclei are discussed. These hRPA and pRPA are formulated 
based on the Hartree-Fock ground state. An extension of hRPA and pRPA based on a correlated 
ground state is given using time-dependent density-matrix theory. Applications to the single-particle 
states around are presented. It is shown that inclusion of ground-state correlation affects ap- 
preciably the results of hRPA and pRPA. The question of the coupling of the center of mass motion 
of the core to the particle (hole) is also discussed. 
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The one-particle states and one-hole states are basic 
excitation modes of a nucleus and other many body sys- 
tems. The experimental data on the properties of the 
single-particle states have been accumulated using nu- 
clear reactions such as one-nucleon transfer, pickup and 
knock-out reactions [1], and it has been found that there 
is a substantial depletion of the spectral strength of the 
single-particle states. Theoretical studies have shown 
that the strong short-range and tensor components of 
the nucleon-nucleon interaction are responsible for a part 
of the depletion [2] and a substantial part of the frag- 
mentation of the single-particle strength is due to the 
coupling to low- lying collective modes [3, 4]. The stan- 
dard approach to study the single-particle properties may 
be the Green's function method. Various theoretical ap- 
proaches have been proposed to implement the coupling 
to low-lying collective modes into the self-energy of the 
Green's function: the particle-phonon coupling model 
[3, 5-7], the Tamm-Dancoff approximation (TDA) [8] and 
the more recent Faddeev random-phase approximation 
(FRPA) [9]. In the present paper we give a formula- 
tion of the hole-state RPA (hRPA) using the equation of 
motion approach (EoM) [10], which has often been used 
to derive the standard RPA, and discuss some aspects 
of hRPA such as the relation to the particle-state RPA 
(pRPA), which have not been clarified so far in the liter- 
ature [8, 11]. We also present an extension of odd A RPA 
(oRPA) based on a correlated ground state obtained from 
the time-dependent density-matrix theory (TDDM) [12- 
14]. The influence of the center-of-mass (c.o.m.) motion 
of the even core on the odd system is also discussed. The 
paper is organized as follows: the formulation of hRPA 
and its extension is given in sect. 2, some properties of 
the extended RPA are also discussed in sect. 2, the results 
obtained for the single-particle states around ^^0 are pre- 
sented in sect. 3, and sect. 4 is devoted to a discussion 
and conclusion section. 



Let us consider a nucleus consisting of A nuclcons and 
assume that the total Hamiltonian H consists of the ki- 
netic energy term and a two-body interaction. Let us 
assume that JO) is the ground state of the A nucleon sys- 
tem with A even and with energy Eq and \fj,) an exact 
eigenstate of the Hamiltonian for the A — 1 system with 
an eigenvalue E^^ {H\fi) = i?^]/i)). 



A. Equations of motion for transition amplitudes 

In direct reaction theories such as the the distorted 
wave impulse approximation and the distorted wave Born 
approximation the differential cross section for one nu- 
cleon transfer reactions is related to the spectral function 



'\aM{fiK\0)Siuj + E^-Eo), (1) 



where and a J are the annihilation and creation opera- 
tors of a nucleon in a single-particle state a, respectively. 
We consider the equations of motion for the transition 
amplitudes and X^^.^ from the A nucleon system to 
the A — 1 nucleon system. These amplitudes are defined 
by 



(o|4Im), 

where : : implies 

Here, n^ai is the occupation matrix given by 

riaa' = (Ola+flalO). 
From the EoM relation 

{0\[H,al]^{Q\al{E^-H), 



(2) 
(3) 

(4) 

(5) 

(6) 
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we obtain the equation for a;^ 

(0|[F,a+]|M)=a;^(0|a+|M)=a;^<, 



(7) 



where ijJ^ = Eq — E^. The commutator on the left-hand 
side of the above equation includes terms with a^a^a^. 
Therefore, cc^ couples to X^^^.^. In a way analogous to 
that used in deriving Eq. (7), we obtain the equation for 



{0\[H,:a+aJa^ :]|m> =i^i^{0\ ■ a+aja^ : \fi) 



(8) 



On the left-hand side of the above equation there appear 
expectation values of the terms consisting of three cre- 
ation operators and two annihilation operators such as 
(0|aJ^a^^oJ^OA4aA5 1/^), which implies the coupling to a 
higher-level amplitude (0| : a+oja+a^/aa' : To close 
the chain of the coupled equations, we factorize these 
terms using and X^^^.^ as 

(0|aAiaA2aA3'^^4aA6lM) ~ -^'5("'A4A3^AiA2:A6 

+ CxMXsX'^x,), (9) 

where the correlation matrix Capa'/S' is defined by 
Cafia'fS' = (0| : a^iOpiapaa '■ |0) and AS{ ) means that 
the terms in the parentheses are properly antisymmter- 
ized [12]. The obtained coupled equations are written 
as 



Ai A2 As 



(AiA2|i;|aA3)X^^A.:A3=0, (10) 



Ai A2 A3 A4 

+ '5A3,fl("-A4Q;"'7A2 + C7A4A2a) 

+ ^A27('^A3tt"A4/3 + 2^^3A4a/3))2;Ai 

+ '5A3a%Ai^A2/3:A4 + '^A3,3n7A2 ^Aia: A4 
- '5Ai7("A4a^^A2:A3 + '^A3/3^aA2:A4 ) 

+ 2^^X30^X4.0 - <^A3C««A4/3 + <^A3 

A2:7-l 

= 0, (11) 

where the subscript A means that the corresponding ma- 
trix is antisymmetrized and the single-particle states are 
chosen as the eigenstates of the matrix 



{a\t\a')+Y,{aX\v\a'\')Anx'x- 



(12) 



AA' 



Here t is the kinetic energy operator. Equations (10) and 
(11) are written in matrix form: 



a c \ X 



(13) 



The matrix elements of the above equation are given 
in Appendix A. The normalization of the amplitudes is 
given by 



Xi" 



(14) 



where x'^*, and X^p*^,^, are the left eigenvector of Eq. 
(13). The occupation matrix and the correlation matrix, 
which enter Eq. (13) and which describe the ground-state 
correlations in the A nucleon system, can be determined 
in the framework of Time Dependent Density Matrix 
(TDDM) theory: the TDDM equations [12, 14] consist 
of the coupled equations motion for riaa' and Ca^a'P', 

i?inaa' = (0|[a+,aa,i?]0), (15) 



ihCaPa'P' 



a^.atapaa :,-H"]|0). 



(16) 



The right-hand side of Eq. (16) contains the expectation 
values of three-body operators, which are approximated 
by the products of Haa' and Capa'p' to close the coupled 
chain of the equations of motion. The ground state in 
TDDM is given as a stationary solution of the TDDM 
equations which satisfies fiaa' = and Ca(}a'[j' = 0. The 
stationary solution can be obtained using the gradient 
method [15]. This method will be used in our numerical 
application given later. 

In the Hartree-Fock approximation (HP), Uaa' — 6aa' 
for hole states and n^a' = for particle states, and 
Ca0a'0' = 0. Keeping in Eq. (13) only the amplitudes 
x'^ and Xf^f^,.^, where p and h refer to a particle state 
and a hole state, respectively, corresponds to the TDA 
equation of odd particle systems [16]. However, even 
within the HF-ground state, Eq. (13) can have all com- 
ponents of A:^^.^; A,^,j,.p, A^p,.^, A^^,.^„, X^,.p„, A^p.p, 
and A'^p.^,. Such equations have been proposed for the 
first time in Ref. [11] and they have been applied in Ref. 
[8]. This very much extended configuration space actu- 
ally leads to some difficulties which have been discussed 
in Ref. [8]. We will take up this discussion again below. 

In the following wc discuss the relation of with naa. 
Using Eq. (10) for x'^, and the complex conjugate of Eq. 
(10) for x^, we can eliminate obtaining an equation 
for E^<'(^S)* 

(e„-e„OE<'«)* 

+ ^ [{aX3\v\X,X,}^x^„,iX^,^,^.,,J* 

A1A2A3 

-{X,X2\v\a'Xs)J2Kxr.Xsi<r] = 0- (17) 
On the other hand the stationary condition fiaa' 

= for 

Eq. (15) gives [15] 

1 A2a'A3 

Ai A2 A3 

-(AiA2|t;|a'A3)C„A3AiA2] = 0. (18) 



Equations (17) and (18) suggest that Y^t^^a'^^a)* and 
^a'/3':/3(^a)* Correspond to riaa' and Capa'd'-, respec- 
tively, though the symmetry of Ca^a'p' under the ex- 
change of a and /3 is lost in 'Yl^j.^a'P'-pi^a)* ■ ^'■^ ^"^^^ 
show below that Uaa = Ylin ^a(^a)* approximately holds 
in the applications to ^^0. 



B. Equation of motion approach with excitation 
operator 

Equation (13) lacks some effects such as self-energy 
contributions in the configuration X'^i^.^, which should 
be included when a correlated ground state is used. In 
order to take account of such effects, we present another 
formulation which is based on EoM [10] . Introducing the 
excitation operator g+ 

a a/37 

and assuming, as usual, (?^|0) = |^) and (j^|0) ~ (for 
the existence of such a relation, see below), we obtain 
from Eqs. (7) and (8) 



A C 



B D I \ V j \ N21 N22 

where the matrices are defined as 



(20) 



A{a:a') = {0\{[H,at],ao^,m, (21) 
B(a/37 : a') = {0\{[H, : a+a+a^ :], a^^lO), (22) 
C{a : a'/3Y) = (0|{[^?,a+ :], : a+a^,a„'}|0), 
D{aj3j : a'13'j' ) 

= {0\{[H , : a+a+a^ :],: a+a0,aa' :}|0), (23) 

iVii(a : a') = (G|{a+, a„,}|0) = 6^^,, (24) 
7Vi2(a : a'/3'j) = (0|{a+, : a+ap,a^, :}|0) = 0, (25) 
N2i{ap-f : a') = (0|{: a+a+a^ :,aa-}|0) = 0, 
N22{al3j : a'l3'j' ) 
= (0|{: a+aja^ , : a+,a^'aa' :}|0). (26) 

Here { } implies the anticommutator, {A, B} = AB + 
BA. The norm matrix N22 is given in Appendix A. The 
matrix elements in Eq. (20) can be expressed using those 
in Eq. (13) such as 



A^ ax Nil 
B = bx Nil, 
C = cx N22. 



(27) 
(28) 
(29) 



The matrix D consists of the two types of terms, one 
expressed by Di = x N22 and the other given by D2, 
which originates from the terms with : a^^a^^a^^a\^a\^ : 




(a) 



(b) 



FIG. 1. (a) Self-energy contribution to a particle state and 
(b) that to a hole state. The ellipses denote C^p^'p' the 
dots the residual interaction. 



in [if, : a^d^a-y :]: 

[ H, : a+o^a^ :] = ^c{al3^ : X)a+, 

a' 

+ E d{af3^ : a'/S'j') : a^,a^,a^> : 

a'0'j' 

+ E e(a/37 : A1A2A3A4A5) 



Ai A2 A3 A4 A5 



(30) 



Using 



A^32 ( A1A2A3A4A5 : al3j) 

= (0|{: «Ai"A2«A3«A5aA4 : a^a^aa :}|0), (31) 

D2 can be expressed as e x A^32. These terms include, for 

example, 

--^^aa'SpB' E ('^1'^2|v|7'A3)aC'7A3AiA2) 

Ai A2 A3 

which is a self-energy contribution to the state 7. The 
self-energy contributions are schematically shown in Fig. 
1. 

The normalization of the amplitudes is given by 
The closure relation is written as 

ir • y ■>{ 

where / is the unit matrix. Wc refer to the formulation 
Eq. (20) as the extended odd-RPA (EoRPA). 



^ = /, (33) 



In the following we discuss the relation between Eqs. 
(13) and (20). The transition amplitudes and X'^^.^ 
are given by and F^^.^ as 

{t')-{ZZW)- 

Inserting this expression into Eq. (13), we obtain 

The difference between Eqs. (20) and (35) and thus be- 
tween Eqs. (20) and (13) resides in the matrix D. Some 
effects of the ground-state correlations such as the self- 
energy contributions arc missing in Eq. (35) and thus in 
Eq. (13) as mentioned above. The importance of these 
missing terms will be discussed below in the application 
section. 



1. Symmetry properties 

First we show that the Hamiltonian matrix of Eq. (20) 
is hermitian. We use the operator identity 

(0|{[ff,i],B}|0) + (0|{[F,B],i}|0) 

= (0|[F,{i,B}]|0). (36) 

In Eq. (20), in the matrix A, the operators A and B are 

identified with a+ and a^i, respectively. Since {A, i?} 
is unity, the right-hand side of Eq. (36) vanishes, which 
means (0|{[if, i], S}|0) = -(0|{[ff,B],i}|0> and 

A{a : a'Y = -{<d\{[H,a^\,a+m, 

= (0|{[if,a+],a„}|0) = A(a':a). (37) 

In the case of the matrix B in Eq. (20) A is : a+a^a-y : 

and B is a„', and {A,B} is reduced to a one-body 
operator. Due to the ground-state condition Eq. (15) 
the right-hand side of Eq. (36) vanishes, which means 
(0|{[ff,i],B}|0) = -(0|{[ff,B],i}|0) and 

= {Q\{[H,a+],:a+a0ao, :}|0) 

= C{a' : ap-f). (38) 

Similarly, for the matrix D in Eq. (20) A is : a+a^a^ : 

and B is : a^itt/s'tta', and {A, B} is reduced to at most a 
two-body operator. Due to the ground-state conditions 
Eqs. (15) and (16) the right-hand side of Eq. (36) van- 
ishes, which implies 

£l(a/37 : a'p'-f')* = D{a'p'i : a^-f). (39) 

Therefore, the Hamiltonian matrix in Eq. (20) is hermi- 
tian. In the applications, shown below, we do not take all 
the matrix elements of Uaa' and Cafja'P'i which causes 




(a) (b) 



FIG. 2. (a) Mass operator for a particle state described by 
^pp'-h ^'^'^ (^) ^^^^ ^^"^ * hole state described by Y^^^,,^. The 
circles mean the propagators given by i^^/.^^ ii^)) ^^'^ ^hh'-p 
((b)), and the dots the residual interaction. 

a violation of the hermiticity of Eq. (20), though it will 
turn out to be small. 

Next we discuss the relation between the formulations 
for a hole state and a particle state. We can obtain a 
formulation for a particle state using the excitation op- 
erator 

9m = E + E ^a/3:7 " " " (^0) 

Oi a/37 

Since this operator is the conjugate of the hole-state ex- 
citation operator Eq. (19), it is easily shown that the 
formulation for a particle state is given as 

{bd) ( I'O = "^'^ ( ivJ! nZ ) ( J'O '^^^^ 

where the superscript t means the transposition of the 
corresponding matrix and is defined by ~ E^^ — Eq. 
Equation (41) imphes that {z^^, Z^) is the left-hand 
eigenvector of Eq. (20). Thus Eq. (20) gives simulta- 
neously the particle states and the hole states. This is 
completely analogous to pp(hh)RPA (see Ref. [16].) 

2. Hartree-Fock approximation for the ground state 

If we make the usual approximation to take for the 
ground state |0) the HF one, N22 in Eq. (26) becomes 

X K7 + <a"°/3 - <7"°" - <74/3)'(42) 

where is equal to 1 or 0. In HF, N22 is non- vanishing 
only for F^'^,.^ and y^'^;,,,^. These amplitudes Y^^.^^ and 

'^hh'-p correspond to the backward amplitudes of yf'^ and 
yjj^, respectively. Hereafter we refer to this formulation 



5 




of N22 are nonvanishing for all configurations due to the 
ground state correlations and, therefore, all other Yj^^.^ 
could be included, in principle. We investigate the effect 
of inclusion of those other amplitudes in EoRPA in some 
limited cases. 




(a) 



(b) 



FIG. 3. (a) Mass operator for a particle state described by 
^hh'-p (^) ^^^^ ^'^^ ^ hole state described by Y^^,,^^. The 
circles mean the propagators given by Yj^^.^ ((a)) and Y^^,,^ 
((b)), and the dots the residual interaction. 



consisting of the four amplitudes, y^, j/^, Y^^.^^ and Yf^^.^ 
as odd-RPA (oRPA). The mass operators of the one- 
body Green's function derived from oRPA are schemati- 
cally shown in Figs. 2 and 3. Since oRPA describes the 
hole states and particle states simultaneously, the single- 
particle strength can be spread over both positive and 
negative energy regions. We consider that the strength 
below the Fermi energy of the core nucleus belongs 
to the states in the A — I system, while that above ep 
to the states in the A -|- 1 system. The TDA hole-state 
equation is obtained by keeping only 7/^ and y^/.^: since 
the coupling of to Y/^f^i.^ is included in addition to the 
coupling to the backward amplitude Y^^i.f^, our oRPA ac- 
tually corresponds to some sort of second RPA for even 
nucleon systems. One may think that in addition the am- 
plitudes y^.^/ and Y^^,„ should be included in oRPA 
because they respectively express the backward propa- 
gations of the particle - hole pair and the hole-hole pair 
in ij^/.p. However, these amplitudes cannot be included 
because the norm of these amplitudes is not defined in 
HF (the matrix elements N22 for Yj^^.^^, and Y^^i.^,, van- 
ish in HF). As mentioned above, the formulation Eq. 
(13) allows us to implement all the X^^.^ amplitudes 



including Xj^^.j^, and X^pi.^,, because there is no restric- 
tion of the norm matrix. (If Eq. (35) is used instead 
of Eq. (13), X'j^^.^i and X'^^,.^,, are projected out, how- 
ever.) The inclusion of all the amplitudes of X'^^.^ can 
give quite unphysical results because the sum of some 
unperturbed energies corresponding to X'^^ fall near 
ei?, which makes it difficult to distinguish the hole states 
from the particle states (see also Ref. [8]). For these rea- 
sons we mainly present the results in EoRPA calculated 
using only the four amplitudes corresponding to y^, y^, 
^pp'-h ^^'^ ^hh'-p oRPA, although the matrix elements 



3. The RPA ground state wave function 

The choice of the subspace spanned by the afore men- 
tioned four amplitudes y'^,y^, Y^^,.^^, Yj^^,.^ may be given 
a different rationale. We consider the following two quasi- 
particle operators which consist only of the forward and 
backward amplitudes 

<it = Y.y>t ^hh':p«' (43) 

p hh'p 

h pp'h 

and neglect the coupling of to ^^^/.^ and that of 
to i^/f^/.p. This oRPA scheme actually corresponds to 
the standard RPA for even nucleon systems. This can 
for example be seen in the following way. It can easily 
be shown that the operators and qp kill the following 
RPA vacuum, i.e. q\Z) = with 

under the conditions 



^^Vp* ^PP'hh' — Yf^^i.p, 



E 



P* TAP* 



(46) 



(47) 



where |HF) is the HF ground state of an even A system. 

These two quasiparticles (one for the particle addition 
(a) and one for the particle removal (p)) span, as seen, ex- 
actly the space of the four amplitudes discussed in II. B. 2. 
However, the single equation for the four amplitudes is 
now split into two independent 2x2 equations corre- 
sponding to the two operators introduced in Eqs. (43) 
and (44), respectively. Using in these equations the HF 
ground state as in II. B. 2, we see, that we have one type 
of 'forward' going amplitudes and one type of 'backward 
going' amplitudes in analogy with what we know from 
standard ph RPA for even systems with corresponding 
amplitudes X and Y. As a matter of fact, it recently 
has been shown [17] that also for the standard ph-RPA a 
generalized operator can be found which annihilates the 
state Eq. (45). It is given by the following form 



Qu - ^[Xphdtap - Yp,^a:l;ah] 

ph 

+ ^ E ^PiP^phO-p^ap^a+ah 

phpip2 

12 ^h,h2phah,ah2apah- 



(48) 



phh^ho, 
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This destruction operator kills the vacuum Eq. (45), i.e. 
Q\Z) =0, under the conditions 



Zphp'h' — ^YphXp'h' 

''ipip2ph — 2 y~!"^pift,i^pp2^fei 

hi 

'nhih2ph ~ n -^pi hi Zppi hh2 ■ 



(49) 
(50) 
(51) 



pi 



We see that there are additional terms to the standard 
ph-RPA operator which contain specific two-body terms. 
The corresponding terms in can schematically be ob- 
tained in augmenting the addition operator Eq. (43) by a 
destructor an and the removal operator Eq. (44) by a cre- 
ator Up . The ?7-terms are also small amplitude (backward 
going) terms which can be added to the standard RPA, 
evaluated with the HF state. They improve the results 
of standard RPA [18]. We, therefore, see that complete 
consistency between RPA in even and odd systems can 
be achieved. 



4- Green's Function Description 

It may be instructive to cast the above amplitude equa- 
tions into Green's function language. For this we write 
down a Dyson equation 

Gkk' = Gl^kk' + ^ M^fc^ fe, , (52) 



where 



1-n 



(0) 



,(0) 



— + 



(53) 



LJ - Sk + ir] CO - Sk -ir] 

is the free or HF Green's function with the occupation 
numbers n'j^^ equal to or 1. The mass operator is given 

by 

Mkk' = E {kh\v\p^P2) 

ahh'pip2p'iP2 

X {p'^p'Mh'k') 



' w - + ir] 



UJ-ila -I 

ppp ri\ii2 ti l '''2 

X {h^h'Mv'k'), 



IT] 



(54) 



where Y"'P and f2„,p are the TDA 2p-lh and 2h-lp 
amplitudes and eigenvalues, respectively, obtained from 
the corresponding TDA equations [16]. In the case 
where we strictly work with oRPA corresponding to 
the ground state Eq. (45), the coupled system of 
particle and hole propagation in above Dyson equation 
decouples into two separate Dyson equations, one for 
the particles (with the a part of the mass operator 
corresponding to Fig. 3(a)) and one for the holes (with 



the p part of the mass operator corresponding to Fig. 
3(b)). It may certainly be appealing to work with an 
approach which is based on a ground state wave function. 

Going beyond the use of a HF ground state, we can do 
as in this work considering EoRPA as described above. 
However, there is also the possibility to mix even and odd 
RPA's. For example it has turned out that Self Consis- 
tent RPA (SCRPA) based on the vacuum Eq. (45) gives 
very good results [17]. For instance, it also solves the 
two particle case exactly. One thus could use SCRPA 
to calculate the correlation functions appearing in an ex- 
tended oRPA. To use the ansatze Eqs. (43) and (44) 
directly seems difficult, since they correspond to a non- 
linear transformation among the fermion operators. 



5. Spurious modes 

First we discuss an RPA-like formulation that can 
bring the c.o.m motion of an odd system at zero excita- 
tion energy. We consider for an A -|- 1 system the ground 
state |<I>o) and an excited state |$^) with excitation en- 
ergy w^. Using the equation of motion 



($o|[aiaa,i?]|<E>^) =a;^($o|aiaa|*^) 



(55) 



and assuming |$o) — 0+ |HF), we obtain the following 
equation 

^I^Ka' = (fa - ea'Xa' 

-«a-<'a')E("^'l^l"'^)^^AA' 



AA' 



MJa^AA' 



+ Y[Sap{p^'\v\a'^)A - 6a'p{aX'\v\pX), 

AA' 

+ Yl^{apM)^P)Ax''^„, - {Xp\v\a'p)Ax''„^]. (56) 

A 

The first two lines of the above equation have the same 
form as the standard RPA for an even A system, and the 
third and fourth terms are due to the additional nucleon 
in a particle state p. For the total momentum operator 



P = Y{a'\ - ihV\a)a+aa, 



(57) 



which satisfies [P,H] = 0, we evaluate w^(i>o|P|^'^) as 



aoi' 



(a'|-^^V|aK<„,. (58) 



Using the right-hand side of Eq. (56) and the trans- 
lational invariance of the interaction [19], we can show 
ujfi{^o\P\^fi) = 0, which implies = 0. Thus the exci- 
tation energy of the c.o.m motion of an odd system given 
by Eq. (56) is zero from the ground state |$o) and Cp 
from |HF). In order to obtain this conclusion, however, 
we need to include all components of x'^^, because of the 
last two terms on the right-hand side of Eq. (56). 
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Now we discuss the e.o.m of a core nucleus in odd A 
nuclei whose treatment is of particular relevance. In the 
standard particle vibration coupling model [3, 5] the spu- 
rious mode is simply discarded, first, for the translational 
mode, on physical grounds but also because the RPA 
amplitudes of a zero mode cannot be normalized. On 
the other hand, e.g. in the case of rotations, it would 
be very important to find a way to include the rota- 
tional mode, since it is a physical state. In order to 
learn something about the coupling of single-particle mo- 
tion and recoil of the core nucleus, we first show that 
a;^(/i|Poa|0) = ea{^i\aa\0) holds in the mean-field ap- 
proximation. Using the complex conjugate of Eq. (8), 
we evaluate oj^{iJ,\Paa\0) such that 

w^(/i|Paa|0) = {n\[Paa,H]\0) 

= (M|[P,ff]aa|0) + (M|P[aa,-ff]|0) 

= {fi\P[a^,H]\0), (59) 

where we use [P,H] = 0. If we use the mean- field 
approximation for [aa,H], that is, [aa,H] = eaUa, then 
we obtain aj^(/i|Paa|0) = ea{fJ-\aa\0), which means 
that the strength \{iJ,\Paa\0)\^ is concentrated at the 
state with = In the general case the mean- field 
approximation is not valid as Eq. (10) indicates. In the 
realistic applications of our oRPA or EoRPA approaches 
shown below, we will, therefore, see that a large portion 
of the strength is distributed to an energy region lower 
than Ca, which can be interpreted as a recoil effect of 
the core nucleus. 

In the past, the question of the spurious modes 
appeared essentially in the particle-vibration coupling 
model [5] which is derived from the Green's function 
method factorizing in the mass operator the 2p-lh (2h- 
Ip) propagator into an ph-RPA propagator and a HP 
single particle propagator. In the spectral representation 
of the RPA propagator the spurious mode is then dis- 
carded because of the zero energy mode and the ensuing 
diverging amplitudes. On the other hand, if one could 
solve the 2p-lh (2h-lp) propagator in the mass operator 
exactly (e.g. in a model) or with a consistent higher order 
theory, surely no problem with a spurious motion of the 
core nucleus would be present. From our analysis above, 
it appears that the mass operator should be calculated 
with 2h-lp (2p-lh) TDA amplitudes. It could very well 
be that this approach gives more realistic results than 
the particle vibration coupling model where the spurious 
mode is discarded. That is what our derivation seems to 
indicate. 

In any case, e.g. in the case of rotations, it would be 
necessary to include this mode, since it is physical. One 
could push the argument even further and assume that, 
since, e.g. the rotation is very collective, the factorization 
of the 2h-lp (2p-lh) TDA into a ph-TDA + plus a hole 
(particle) is a good approximation (the neglected terms 
coming only from exchange). Because of its strong col- 
lectivity, eventually all the other couplings to intrinsic ph 
modes could be neglected. Actually analogous questions 



would arise in cold fermionic atom systems where one 
could ask the question what happens to an odd fermion 
which is coupled to the so-called Kohn mode, i.e. a coher- 
ent c.o.m. motion of the underlying even system, in the 
external harmonic container. Since the mass of the core 
can be very large, e.g. with a million of atoms, the fac- 
torization can become quite valid and also the ph-TDA 
for the Kohn mode will become very collective. It could 
be interesting to investigate this question in more detail 
theoretically and experimentally because the treatment 
of Goldstone modes in single-particle mass operators is, 
to the best of our knowledge, an unsolved problem. 



C. The A = 2 case 

We show that our formulation is exact for an A = 2 
system. In the case of an A = 2 system, TDDM gives the 
coupled equations of motion for n^a' and the two-body 
density matrix pafsa' p' , which are defined as 

n^c.'{t) = {Ht)\al'Ci,,\^{t)). (60) 
Pai}a'/}'{t) = {^{t)\a+,a-^,ai3aa\^{t)}, (61) 

where |<I'(f)) is the time-dependent total wavefunction 
\^{t)} = exp{-iHt]\<!>{t = 0)). The equations in TDDM 
are written as [20] 

ihhaa' = ^{{a\t\X)nxa' - {X\t\a')nax) 

Ai A2 A3 

- PaAiA2A3(A2A3|?^|a'Ai)], (62) 

ihpal3a'l3' = ^((a|t| A)/9A/3a'/3' + {P\t\>) PaXa' 13' 
A 

- {X\t\a')pal3\l3' - {X\t\P')poi0Oi'x) 

+ ^[i.al3\v\\i\2)pX-,X^oc'0' 

X1X2 

- (AiA2|w|a'/3')Pa/3AiA2]- (63) 

Here the single-particle states are arbitrary. Since there 
are no higher-level reduced density matrices in an A = 2 
system, these two equations are exact. When the two- 
body density matrix in Eq. (62) is approximated by anti- 
symmetrized products of the occupation matrices. Eq. 
(62) is equivalent to the equation in the time-dependent 
HF theory. The ground state is given as a stationary 
solution of these equations. 

The equation for the transition amplitude is 

^m\a)-6^xuJ,)x';i+ Yl (AiA2b|aA3)Xj:^A2:A3 

A A1A2A3 

= 0, (64) 
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where Xj^^,^ = (0|a+a^a^|/x). The equation for X^^.^ is 
given as 

A 

Ai A2 

= 0. (65) 

Since there are no higher-level transition amplitudes in 
an A = 2 system, these two equations are also exact. 
From Eq. (64) we obtain 

fXOLCX.' 

= ^{a'\t\a)naa' 

OlOl' 

+ 2 X] ('^1'^2|i'|A3A4)PA3A4AiA2 

Ai A2 A3 A4 

= (0|H|0), (66) 
where Uaa' and are exactly given by 

Pa0a'0'=J2K'0':0i<r- (68) 

Equation (66) corresponds to the relation between the 
total ground-state energy and the single-particle Green's 
function [21]. 



TABLE I. Single-particle energies and occupation prober- 
bilities riaa calculated in TDDM. The single-particle energies 
in HF are given in parentheses. 





€c [MeV] 


Ticket 


orbit 


proton neutron 


proton neutron 


lSl/2 
lP3/2 
lPl/2 
1^5/2 
2Sl/2 


-32.5 (-32.1) -36.2 (-35.9) 
-18.3 (-18.2) -21.8 (-21.8) 
-12.3 (-12.0) -15.8 (-15.6) 
-3.8 (-3.8) -7.1 (-7.2) 
1.2 (1.5) -1.5 (-1.2) 


0.98 0.98 
0.93 0.93 
0.91 0.91 
0.08 0.08 
0.02 0.02 



B. Ground state 

The occupation probabilities calculated in TDDM are 
shown in Table I. The largest deviation from the HF 
values (ti^„=1 or 0) is about 10%, which means that the 
ground state of -"^^O is a strongly correlated state. A re- 
cent shell-model calculation by Utsuno and Chiba [22] 
also gives a similar result for the ground state of ^^O. 
The correlation energy Ec in the ground state, which is 
defined by Ec = I]a/3a'/3'("/3|wla'/3')Ca'/3'a/3/2, is -19.6 
MeV. A large portion of the correlation energy is com- 
pensated by the increase in the mean-field energy due 
to the fractional occupation of the single-particle states. 
The resulting energy gain due to the ground-state corre- 
lations, which is given by the total energy difference be- 
tween HF and TDDM, is with 5.2 MeV relatively small. 
Such kind of scenario is similar to the one well known 
from BCS theory [16]. 



III. APPLICATIONS TO ^''O 

A. Calculational details 

In this paper, we make a first schematic application of 
our theory to proton hole states in ^^N and proton par- 
ticle states in ^^F. We do not consider the correspond- 
ing neutron states because there are less experimental 
data. We consider the lsi/2, lp3/2, lpi/2, lc^5/2, 2si/2, 
1^3/2, 2^3/2, 2pi/2 1/7/2 and I/5/2 states for both pro- 
tons and neutrons. The continuum states are discretized 
by confining the single-particle wavefunetions in a sphere 
of radius 12 fm. We use a simplified residual interaction 
which consists only of the to and terms of the Skyrme 
HI force. Its strength is reduced by 20% to put the spu- 
rious c.o.m motion of at approximately zero energy 
in the standard RPA. For the ground-state calculation 
of ^^O in TDDM, we only use the bound single-particle 
states, the lsi/2, IP3/2! lPi/2! lc^5/2 and 2si/2 states, and 
consider only the two particle - two hole type correlations 
in Ca/3a'^'- We also neglect the off-diagonal elements of 
Uaa' between the lsi/2 and 2si/2 states. The am- 
plitudes for the proton 2si/2, 2pi/2 and 2pi/2 states are 
neglected because their contributions are negligible. 



C. Spectral functions 

In Figs. 4 and 5 the spectral functions of the proton 
lpi/2 and lp3/2 hole states in ^^N calculated in EoRPA 
(Eq. (20)) (solid line) are shown, respectively, and com- 
pared with the results in TDA (dotted line). Since the 
results in oRPA are similar to the EoRPA resiilts, they 
are not shown. As already mentioned, in the EoRPA cal- 
culations we consider only the same Y^^^.^ amplitudes as 
those used in oRPA. To facilitate a comparison of various 
calculations, we smooth the distributions using an arti- 
ficial width r = 0.5 MeV. As shown in Table I the HF 
energies of the proton lpi/2 and IP3/2 states are —12.0 
MeV and -18.2 MeV, respectively. In TDA the main 
peak is shifted upwards from the HF position due to the 
coupling to the configurations Y/^f^,.^ whose unperturbed 
energies are distributed below —40 MeV. In EoRPA (and 
oRPA) the main peak is slightly shifted downwards from 
the HF position due to the coupling to both V^/.^ and 
the backward amplitudes ^^/.^ whose unperturbed en- 
ergies are located above MeV (see Fig. 4 and Fig. 
5). The strength distribution in the positive energy re- 
gion corresponds to the states in -'^^F. The strengths of 
the main peak of the proton lpi/2 state calculated in 
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FIG. 4. Spectral function of the proton lpi/2 state in ^""N 
calculated in EoRPA (solid line). The dotted line shows the 
result in TDA. The distributions are smoothed with an ar- 
tificial width r = 0.5 MeV. The strength distribution in the 
positive energy region is due to the coupling to the backward 
amplitude Yj^,. and indicates the states in ^'^F. The small 
strength distributions in the positive and negative energy re- 
gions are shown in the insets. 
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FIG. 5. Same as Fig. 4 but for the proton lp3/2 state. 



EoRPA, oRPA and TDA arc 0.88, 0.82 and 0.95, respec- 
tively. When the forward amplitude Yj^y is neglected 
in oRPA, the main peak is further shifted cfown to —14.9 
MeV and has strength 0.89. This indicates that the cou- 
pling to the backward amplitudes i^^/.^j plays an impor- 
tant role in depleting the single-particle strength. The ef- 
fects of the ground-state correlations included in EoRPA 
play a role in slightly reducing the correlations in oRPA 
due to fractional occupation of the single-particle states. 
The sum of the strength of the proton lpi/2 state dis- 
tributed in the negative energy region is 0.91 in EoRPA, 
which corresponds to naa = 0.91 in TDDM, (sec Table 



FIG. 6. Same as Fig. 4 but for the proton lsi/2 state. 



good approximation. The single-particle strengths of the 
main peak of the proton 1^3/2 state calculated in EoRPA, 
RPA and TDA are 0.88, 0.82 and 0.93, respectively. The 
sum of the occupation probabilities of the proton IP3/2 
state distributed in the negative energy region is 0.93 
in EoRPA, which corresponds to Haa = 0.93 in TDDM. 
Summing the whole spectral weights in negative and pos- 
itive energy regions gives, of course, the sum rule value 
of one. 

The results for the proton lsi/2 state are shown in 
Fig. 6. The HP energy of the proton lsi/2 hole state 
is —32.1 MeV. The strength is fragmented due to the 
coupling to the configurations Y/^f^,,^: the unperturbed 

energy of the configuration (lpi/2)~^(lp3/2)~"'^l'^5/2 is 
about -26 MeV. Since the backward configurations Y^,,f^ 
are energetically well separated, there is no significant 
difference between the TDA and oRPA results. There- 
fore, the oRPA result is not shown in Fig. 6. Comparing 
with the results obtained from Eq. (35), we found that 
the D2 term in the matrix D, which is given by e x iV32 
and describes the self-energy contributions to the one 
particle - two hole configurations, play a role in shifting 
the strength to lower energy region. The summed occu- 
pation probability of the proton lsi/2 state distributed 
in the negative energy region is 0.98 in EoRPA, which 
corresponds to n„Q = 0.98 in TDDM. 

The spectral function of the proton 1^5/2 state in ^^F is 
shown in Fig. 7. The HF energy of the 1^5/2 state is —3.8 
MeV. The main peak is shifted downwards from the HF 
position in TDA due to the coupHng to i^^/./j, while, on 
the contrary, it is shifted upward in oRPA due to the ad- 
ditional coupling to the backward amplitudes Yj^^,,^. The 
ground-state correlations included in EoRPA play a role 
in slightly reducing correlations in oRPA. The states lo- 
cated below the single-particle energy of the proton lpi/2 
state correspond to the states in ^^N. The summed occu- 
pation probability of the proton 1^5/2 state distributed 



I). Thus the relation n^a = X^u K^I'^alO)!^ holds to a below the proton lpi/2 state is 0.06 in EoRPA, while the 
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FIG. 7. Same as Fig. 4 but for the proton 1^5/2 state in ^'^F. 
The dot-dashed hne shows the result in oRPA. The strength 
distribution below —15 MeV is due to the coupling to the 
backward amplitude Yt^^,. and shows the states in ^^N. 



FIG. 9. Spectroscopic factors for the proton lsi/2 and 2si/2 
states calculated in oRPA are compared with experiment [23] 
(red bars). 
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FIG. 8. Spectroscopic factors for the proton lpi/2 and IP3/2 
states calculated in EoRPA are compared with experiment 
[23] (red bars). 

corresponding value for Uaa in TDDIVI is 0.08. 

D. Comparison with experiment 

The spectroscopic factors (defined by (2j -t- 
1) xtransition strength) calculated in EoRPA for 
the proton lpi/2 and IP3/2 states are compared with 
experiment [23] (red bars) in Fig. 8. The main peali 
of the proton lpi/2 state is considered as the ground 
state of ^^N and the hole-state energy is measured 
from this threshold in the following. The results in 
EoRPA are reasonable though they overestimate the 
experimental data and cannot reproduce the strength 
distribution around —10 ]\/[eV. This is a common feature 



FIG. 10. Same as Fig. 9 but for EoRPA. 



of TDA and RPA-type calculations [8, 9]. Studies on 
the effect of short-range correlations have predicted a 
strength reduction of about 10% in ^^O [24-26]. The 
spectroscopic factors for the proton lsi/2 and 2si/2 
states calculated in oRPA is compared with experiment 
[23] (red bars) in Fig. 9. The EoRPA results are 
also compared with experiment (red bars) in Fig. 10. 
Since the inclusion of ground-state correlations causes a 
downward shift of the strength, the agreement with the 
data becomes somewhat worse in EoRPA. The strong 
fragmentation below —15 MeV cannot be reproduced in 
these oRPA and EoRPA calculations. Probably higher 
configurations are needed. The spectroscopic factors for 
the proton 1^5/2 and 1^3/2 states calculated in oRPA 
are compared with experiment [23] (red bars) in Fig. 11. 
The EoRPA results are also compared with experiment 
in Fig. 12. Due to a downward shift of the strength, 
the agreement with the data is worsened in EoRPA. We 
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FIG. 11. Spectroscopic factors for the proton 1^5/2 and ld3/2 
states calculated in oRPA are compared with experiment [23] 
(red bars). 
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FIG. 12. Same as Fig. 11 but for EoRPA. 



point out that there is a similar situation in the first 3^ 
state in -'^^O. The effects of the ground-state correlations 
can be included into the standard RPA using riaa' and 
Cafsa'p' as in SCRPA [27]. The first 3" state of ^^^O 
calculated in this modified RPA scheme comes about 
5 MeV higher than the result in the standard RPA. 
This is the same situation as the EoRPA results shown 
above. We found that inclusion of the coupling of the 
particle-hole amplitude to higher two-particle two-hole 
amplitudes brings down the first 3^ state to the right 
position [28]. Therefore, more elaborate calculations 
using a larger number of the Y^^.^ amplitudes and 

also the higher amplitudes Yf^a^,.\M could shift the 



the higher amplitudes Y^^^.^y 
strength upward and bring a better agreement with 
the data. Globally, one may say that the agreement of 
spectroscopic factors with experiment is only marginally 
satisfactory indicating the need for inclusion of higher 
configurations. 



E. Effects of other amplitudes 



We investigate the effects of inclusion of the amplitude 
^hp-h' i'^ EoRPA, which describes backward scattering of 



hp 

a particle - hole pair in Yj^j^i.^. 



We use for Yj^ 



the 



same truncated single-particle space as that used in the 
ground-state calculation since it is important to include 
the self-energy contribution to all single-particle states in 
^hp-h'- reduce the dimension size, we neglect the am- 
plitude The obtained result for the proton lpi/2 
state is shown in Fig. 13 and compared with the result 
(red bars) of the calculation based on Eq. (13) where the 
ground state is assumed to be the HP ground state and 
only the amplitudes X'^p.}^, and Xf^j^,.^ are included. As 
shown in Fig. 13, the inclusion of Xf^^.j^, gives quite un- 
physical results : the main peak is fragmented and some 
states have negative strength. The reason for the frag- 
mentation of the main peak is that unperturbed energies 
of some Xj^j^.f^, fall near the energy of the proton lpi/2 
state. For example, the unperturbed energy of the config- 
uration (lsi/2)^^2si/2 (1^3/2)^^ that couples to the pro- 
ton lpi/2 state is —12.4 MeV, which is close to the energy 
of this state (eq = —12.0 MeV). These unphysical prop- 
erties are not seen in the EoRPA result. We consider that 
this is due both to the self-energy insertion to the config- 
urations Y/^p.f^, and to their small normalization iV22. The 
energy of the configuration F^.^/ is significantly shifted 
by the amount determined by the self energy and the 
normalization. This shift probably plays a role in reduc- 
ing the coupling to the single-hole state. We performed a 
similar EoRPA calculation for the proton lsi/2 state, see 
Fig. 14 and the obtained result (solid line) is compared 
with the EoRPA result without Y/^^.^^, (dotted hne). The 
coupling to Yj^p.j^i plays a role in shifting some strength 
upward, which improves the agreement with the exper- 
iment. However, we found that the inclusion of other 
amplitudes such as Yf^f^,.f^„ and Y/^^.^, brings unphysical 
fragmentation of the strength of the lpi/2 state as seen 
in Fig. 13. Therefore, it requires further investigation 
whether the amplitudes of Y^^^.^ with small normaliza- 
tions should be included or not in EoRPA. 



F. Center of mass motion of ^"O 

Finally we discuss the coupling of a hole state to the 
c.o.m. motion of the core nucleus ^^O using oRPA. The 
strength distribution of |(/j,|PzaQ|0)p among the states 
which couple to the proton lsi/2 state is shown in Fig. 
15, where the proton lpi/2 and 1^3/2 states are taken for 
a. The upper part of Fig. 15 shows the strength of the 
proton lsi/2 state (the same as Fig. 9) and the lower 
part — |(/z|P2aQ|0)p for a = lpi/2 (red bars) and IP3/2 
(green bars). The strength |(/i|Pzaa|0)p for a = lpi/2 
is concentrated in a single state and the coupling of the 
proton lsi/2 state to this state is negligible. Therefore, 
this state may be interpreted as a spurious mode consist- 
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FIG. 13. Strength distribution of the proton lpi/2 state calcu- 
lated in EoRPA including the amplitude Y^^.^^,. The red bars 
show the result of the calculation based on Eq. (13) where 
only the amplitudes X^^,^^, and X'^^.^ axe included under the 
assumption of the HP ground state. Small strengths in the 
low-energy region are shown in the inset. 



FIG. 15. Strength distribution of |{/i|Pzaa|0)p in the states 
which couple to the proton lsi/2 state. The upper part shows 
the strength of the proton lsi/2 state and the lower part 
— |(/i|Pzaa|0)p (in arbitrary units) for a — the proton lpi/2 
state (red bars) and the proton IP3/2 state (green bars). The 
strength |(/^|Pzaa|0)p for a = IP3/2 is fragmented due to the 
coupling to the configurations which have the particle-hole 
pairs with angular momentum L — 2 h. 
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FIG. 14. Strength distribution of the proton lsi/2 state cal- 
culated in EoRPA with (solid line) and without (dotted line) 

hp:h' ' 



ing of a pure c.o.m motion of '^^O and the proton lpi/2 
state, though it is located about 7 MeV below the single- 
particle energy of the proton lpi/2 state (—12.0 IMeV). 
This energy shift from the single-particle energy is related 
to the fact that the TDA calculation for the c.o.m motion 
of ^^O gives the excitation energy of 7.7 MeV. This fact 
may look perturbing, since we liuow that the spurious 
mode of the core comes at zero energy in the standard 
RPA. As already mentioned above, it is not guaranteed in 
an odd system that the c.o.m motion of the core comes at 
zero energy, and the coupling of the spurious mode to the 
physical spectrum is very week, so that the position of the 
spurious mode is not so perturbing. We also performed a 



FIG. 16. Strength distribution of |(/.t|Pzaa|0)p in the states 
which couple to the proton lpi/2 state. The upper part shows 
the strength of the proton lpi/2 state and the lower part 
— \{fi\Pzaa\0)f (in arbitrary units) for a — the proton lsi/2 
state (red bars). 



TDA calculation for the c.o.m motion of -^^^Pb using the 
single-particle states and the residual interaction which 
put the spurious mode at zero energy in RPA and found 
that it comes at 4.7 MeV. This suggests that the c.o.m 
of very heavy systems could come close to zero excita- 
tion energy even in TDA. The strength |(/i|P2aa|0)p for 
a = lp3/2 is rather fragmented. Since the dominant com- 
ponents of the c.o.m motion of the core consists of the 
transitions from the IP3/2 to 1^5/2 states, the large frag- 
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FIG. 17. Same as Fig. 16 but for the proton lp3/2 state. 



mentation of K/ilPzaajO)!^ for a = lp3/2 is explained by 
the coupling to the configurations with different coupling 
schemes of angular momenta: in the case of a = lp3/2 5 
'^ah-p consisting of the particle - hole pairs (lp3/2)^^lc^5/2 
with angular momentum L — 1 h can couple not only 
to Y^j^.^ which has the same particle - hole pairs with 
L = 2 h but also to Y^f^.^ consisting of the particle - 
hole pairs (lpi/2)~^l<^5/2 with L = 2 h. In the case of 
a = lPi/2 the particle - hole pairs (lp3/2)~^lrf5/2 in Y^^.^ 
can have only L — 1 h and does not couple to nearby con- 
figurations. The exchange effect may also play a role in 
weakening the coherence of the c.o.m motion of ^^O for 
a — lp3/2- we tried an oRPA calculation for a = IP3/2 
where all exchange terms are neglected and observed the 
appearance of such a coherent state as that observed for 
a — lpi/2- Thus, in the case of the proton lsi/2 state 
and a = IP3/2 the c.o.m motion of the core nucleus is 
embedded in the physical states of the A—\ nucleus and 
cannot be neglected. The distributions of \{^\Pzaa\Q)\'^ 
among the states which couple to the proton lpi/2 and 
lp3/2 states are shown in Figs. 16 and 17, respectively, 
where a is the proton lsi/2 state. In the case of the pro- 
ton lpi/2 state the strength |(/x|PzaQ|0}|^ is concentrated 
in a single state and the coupling of the proton lpi/2 state 
to this state is negligible. As in the case of the proton 
lsi/2 state this state may be interpreted as a spurious 
mode consisting of the c.o.m motion of ^^O and the pro- 
ton lsi/2 state, though it is located about 5 MeV below 
the single-particle energy of the proton lsi/2 state (—32.1 
MeV). The fragmentation of \ {^\Pzaa\Q)\'^ for the proton 
lP3/2 state is larger than that for the proton lpi/2 state. 
This is explained by the coupling to the configurations 
with different angular momentum couplings: in the case 
of the lp3/2 state the particle - hole pairs (lp3/2)~^ld5/2 
in i^/j.p can carry angular momentum L — 1 h and 2 h, 
whereas the pairs cannot have L = 2 ?i in the case of the 
lpi/2 state. 



IV. SUMMARY 

In this paper, we took up the old subject of the RPA 
approach to odd particle systems. Those equations 
based on the usual equation of motion method (EoM) 
encountered in the past some difficulties [8]. This gave 
raise to the so- called Faddeev-RPA (FRPA) approach 
[9]. However, whenever the RPA breaks down, so does 
FRPA. We located some of the difficulties of the old 
odd particle RPA (oRPA) and proposed some cure, 
limiting the configuration space to the normalizable 
subspace. We showed that p-RPA and h-RPA equations 
give identical results what is very similar to the property 
of pp(hh)RPA for even systems [16]. We also discussed 
the influence of the c.o.m. motion of the core on the 
odd particle (p or h). No difficulty with a break down 
seems to arise. It turned out that the recoil of the 
core influences the spectrum. This aspect may be 
most important for rotational states in deformed nuclei 
where the so-called spurious modes are, in fact, physical 
states. We also showed how to include ground state 
correlations explicitly in EoRPA, similar to what is done 
in TDDM, on top of the oRPA equations. We made a 
first schematic application, using a simplified Skyrme 
force, to the hole and particle states around ^^O. We 
compared Tamm Dancoff, oRPA, and EoRPA solutions. 
It was shown that in some cases all three approaches give 
very similar results but that in others the infiuence of 
extra RPA correlations were significant. The comparison 
with experiment is sufficiently encouraging to develop 
this kind of RPA approach further. In fact, the spirit 
of oRPA is quite close to second RPA. We encountered 
problems for the odd systems e.g. that the spectrum 
becomes too much shifted downwards. Such open 
problems may be a subject for the future. Also the 
connection between a common RPA vacuum in the even 
and odd systems, as proposed recently [17], may be an 
interesting further line of research. 



Appendix A 



a{a : a') — taSaa' (^^1) 
b{al3j : a') ~ ''^^{a' X\v\al3) An-yx 

X 

- '^[{a'X'\v\aX)Ainxi3n^x' +C^xx'i3) 

XX' 

+ {a'X'\v\Xl3)A{nxan^x' +C^xx'a) 
-(a'71 V \XX')A{nxanx'i3 + ^Cxx'ap)], (A2) 

c{a:a'l3'j') = {a'f3'\v\a-f') (A3) 
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d ( a/37 : Q-'P'^') = (e^ + - e^)5aa'5p0'5^^i 
+ ^[{>^oi'\v\a'y')Anryx50i3, - (Aa'|u|/37')A%A^a/3' 

A 

+ {'^P'\v\X^')An\a5l3a' - {l S' \v\\"f') AnxjlSaa' 

- ^5^^\{a'p'\v\a\)Anxfi + {a' l3'\v\\l3)Anxa)]- (A4) 
The norm matrix N22 is given as 

+ ('^7'7^/3a' + C^ip^ai) 

+ ^/3a' {n"y'jnafi' + C^'a-yP')- (A5) 
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